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FOREWORD 


Stability  properties  of  a  charged  beam  propagation  through  a  relativistic 
hollow  electron  beam  are  investigated,  in  connection  with  present  experimental 
applications  in  the  collective  particle  accelerator.  The  stability  analysis  is 
carried  out  for  long  axial  wavelength  and  low-frequency  perturbations.  A  closed 
algebraic  dispersion  relation  for  coupled  transverse  oscillations  is  obtained  for 
the  solid  and  hollow  beams  with  sharp-boundary  density  profiles.  One  of  the  most 
important  features  in  the  analysis  is  that  the  typical  growth  rate  of  the 
transverse  oscillation  is  order  of  the  hollow  beam  diocotron  frequency  wn, 
thereby  severely  limiting  the  solid  beam  propagation  through  a  relativistic 
hollow  electron  beam.  However,  for  a  solid  beam  with  a  small  radius,  the 
fundamental  mode  perturbation  (i.e.,  the  dipole  oscillation)  is  the  most 
unstable  mode. 


IRA  M.  BLATSTEIN 
By  direction 
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INTRODUCTION 

In  recent  years, there  have  been  numerous  theoretical  investigations 

1  2 

of  the  equilibrium  and  stability  properties  ’  in  intense  relativistic 

hollow  electron  beams,  motivated  by  a  variety  of  applications, 

including  microwave  generation  and  amplification,  ’  and  collective 

effect  acceleration5  ^  by  relativistic  hollow  electron  beams.  Moreover, 

8  9 

recent  renewed  interest  in  the  equilibrium  and  stability  properties  ’ 

of  charged  particle  beams  propagating  through  a  relativistic  electron 

beam  originates  from  several  diverse  research  areas,  including 

confinement  and  transport^  ’  ^  of  a  nonneutral  electron  beam,  Che 

collective  ion  acceleration1^  ^  and  electron  beam  propagation15,1^ 

through  a  background  plasma.  In  particular,  the  collective  pat.  tide 
5  6 

accelerator  ’  is  the  acceleration  of  a  charged  particle  beam  propagating 
in  a  spatially  modulated  magnetic  field.  In  this  regard,  this  paper 
examines  the  equilibrium  and  stability  properties  of  a  solid  charged 
beam  propagating  through  an  intense  relativistic  hollow  electron 
beam,  in  connection  with  present  experimental  applications  in  the  collec¬ 
tive  particle  accelerator. 

Investigation  of  the  coupled  transverse  oscillation  is  carried 
out  for  an  intense  charged  particle  beam  that  is  infinite  in  axial 
extent  and  propagates  through  a  relativistic  hollow  electron  beam. 

The  stability  analysis  is  calculated  within  the  framework  of  a  hybrid 
(Vlasov-fluid)  model  in  which  the  hollow  beam  electrons  are  described  as 
a  macroscopic,  cold  fluid  immersed  in  an  axial  magnetic  field  B^e^ 
and  the  solid  beam  particles  are  described  by  the  Vlasov  equation. 

We  assume  that  the  radius  of  the  solid  beam  Rs  is  less  than  the 
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inner  radius  of  the  hollow  beam  and  that  the  hollow  electron 

beam  is  tenuous.  Moreover,  it  is  further  assumed  that  v  /y  «  1, 

s  s 

where  the  subscript  denotes  solid  beam  particles,  is  Budker's  parameter, 
2 

and  Ysn»sc  is  the  characteristic  particle  energy  for  the  solid  beam. 

Equilibrium  properties  and  the  basic  assumptions  are  presented  in  Sec.  II. 

Stability  analysis  of  coupled  transverse  oscillations  is  carried  out 
in  Sec.  Ill  for  long  axial  wavelength  and  low-frequency  perturbations. 

A  closed  algebraic  dispersion  relation  [Eq.  (28)]  of  the  coupled 
transverse  oscillation  is  obtained  for  the  solid  beam  described  by 
Eq.  (6)  and  the  hollow  beam  with  density  profile  in  Eq.  (9).  Equation 
(28)  is  one  of  the  main  results  of  this  paper,  and  can  be  used  to 
investigate  stability  properties  for  a  broad  range  of  system  parameters 
of  experimental  interest.  One  of  the  most  important  features  in  this 
analysis  is  that  the  coupling  between  the  hollow  and  solid  beam  modes  is 
minimum  whenever  both  beams  have  a  same  axial  velocity.  As  a 
particular  example,  in  the  absence  of  either  the  solid  beam  or  the 
hollow  beam,  we  recover  the  previous  results1’1^  from  the  dispersion 
relation  in  Eq.  (28). 

Making  use  of  the  dispersion  relation  in  Eq.  (28),  stability 
properties  for  coupled  transverse  oscillations  between  the  solid 
electron  beam  and  the  hollow  electron  beam  is  investigated  in  Sec.  IV. 
Several  points  are  noteworthy  in  the  analysis  in  Sec.  IV.  First,  for 
a  solid  beam  satisfying  Rs  <<  R^,  the  fundamental  mode  perturbation 
is  the  most  unstable  mode.  Second,  the  typical  growth  rate  of  the 
transverse  oscillation  is  a  substantial  fraction  of  the  hollow  beam 
diocotron  frequency  thereby  severely  limiting  the  solid  beam 
propagation  through  the  hollow  electron  beam.  The  growth  rate  of  the 
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conventional  diocotron  instability  is  much  less  than  that  of  the  coupled 
transverse  oscillation  for  counter-streaming  beams.  However,  the 
growth  rate  of  instability  reduces  drastically  with  decreasing  value 
of  the  density  ratio  ^/n^  of  the  solid  beam  to  hollow  beam.  In  this 
regard,  a  solid  beam  with  a  substantially  reduced  density  can  propagate 
through  a  counter-streaming  hollow  beam  without  a  severe  growth  of 
unstable  oscillations. 

The  ion  resonance  stability  properties  of  a  solid  ion  beam 

propagating  through  a  hollow  electron  beam  is  investigated  in  Sec.  V 

for  a  low  density  case.  For  a  case  of  high  density  ion  beam,  we  urge  the 

8  9 

reader  to  review  the  previous  studies.  *  It  is  found  from  the 
analysis  in  Sec.  V  that  the  maximum  growth  rate  of  ion  resonance 
instability  is  many  times  of  the  diocotron  frequency.  Moreover,  the 
instability  growth  rate  increases  significantly  as  the  solid  beam 
radius  Rg  is  increased. 
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THEORETICAL  MODEL 


The  equilibrium  configuration  consists  of  an  intense  charged 
particle  beam  that  is  infinite  in  axial  extent  and  propagates  through  a 

A 

relativistic  hollow  electron  beam  with  axial  velocity  6^  c  ,gz-  The 
axial  velocity  and  radius  of  the  solid  beam  are  denoted  by  8SC,^Z  ®nd  Rs» 
respectively.  The  inner  and  outer  radii  of  the  hollow  beam  is  represented 
by  and  R2,  and  this  hollow  beam  is  located  inside  a  cylindrical 
conducting  wall  with  radius  R^.  For  simplicity,  we  assume  that  the 
radius  of  the  solid  beam  is  less  than  the  inner  radius  of  the  hollow 
beam,  i.e.,  R  <  R^.  The  applied  axial  magnetic  field  B^z  provides 
confinement  of  the  equilibrium  configuration  in  the  radial  direction. 

In  the  theoretical  analysis,  we  introduce  a  cylindrical  polar  coordinate 
system  (r,0,z). 

In  order  to  make  the  analysis  tractable,  the  following  simplifying 
assumptions  are  made; 

(a)  The  motion  of  the  beam  particles  is  predominantly  in  the  axial 
direction,  and  that  the  transverse  momentum  is  small  in  comparison  with 
the  characteristic  axial  momentum,  i.e.. 


2 

Pr 


<< 


(1) 


where  £  ■  (pr,  p^,  pz)  is  the  particle  momentum  in  laboratory  frame. 

(b)  The  theoretical  analysis  is  carried  out  for  a  tenuous 
hollow  electron  beam  satisfying 


-2  _  2 
“Ph  “  wch  * 


(2) 


.2  2. 

where  £1^  ”  ^7re  ^t/Yhm  P^asina  frequency-squared  of  hollow  beam 
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electron,  »  e®0^hmc  is  the  hollow  beam  electron  cyclotron  frequency, 
and  y. 


(i  -  s*r1/2. 


(c)  It  is  further  assumed  that 


v 

-S.  , 

Ys 

where  s  denotes  solid  beam 


,2  2 
Z  e 


s  2  y 

V  3 


~  «  1  » 


(3) 


particles,  v  is  Budker’s  parameter,  N  *  2ir  dr  r  nu(r)  is  the 

s  s  J0  s 

number  of  particles  per  unit  axial  length,  n^(r)  is  the  equilibrium 

particle  density,  -e  is  the  electron  charge,  Z  is  the  charge  state 

2 

of  the  solid  beam  particles,  Ysfflsc  is  the  characteristic  particle 
energy  for  solid  beam  and  c  is  the  speed  of  light  in  vacua . 

The  analysis  is  carried  out  within  the  framework  of  a  hybrid 
(Vlasov-fluid)  model  in  which  the  hollow  beam  electrons  are  described 
as  macroscopic,  cold  fluid  immersed  in  an  axial  magnetic  field  B^|z 
and  the  solid  beam  particles  are  described  by  the  Vlasov  equation. 

In  this  context,  the  equation  of  motion  and  the  continuity 
equation  for  the  hollow  electron  fluid  can  be  expressed  as 


<fF+*  *  V>^--e(lT+ V^)’ 
It  nh  +  z  •  - 0  • 


(4) 

(5) 


where  and  B^  are  the  electric  and  magnetic  fields,  respectively, - 

n^  and  V  are  the  density  and  mean  velocity,  respectively,  of  the  hollow 

electron  fluid  element,  m  is  the  rest  mass  of  electrons,  and  y(£,t)  * 

2  -1/2 

(1  -  V  •  %/c  )  .  For  present  purposes,  we  assume  that  the  equilibrium 

distribution  function  for  the  solid  beam  particle  is  of  the  form  * 
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fs(^’^)=(fls/2irYsms)<5(H  ~  uspe  "  Ysmsc2) 


X  6  (P  -  Y_m  6„c)  , 

z  s  s  s 


(6) 


22  22  1/2 

where  ftg  is  the  particle  density  at  r  -  0,  H  ■  (mgc  +  c  £  )  + 

eg4>o(r)  is  the  total  energy,  Pg  -  r[pg  +  (eg/c) (rB^/2) ]  is  the  canonical 
angular  momentum,  P_  -  p„  +  (e  /c)A^(r)  is  the  axial  canonical 
momentum,  eg  is  the  charge  of  the  solid  beam  particles,  <f>g(r)  is  the 
electrostatic  potential  for  the  equilibrium  self-electric  field, 

A^(r)  is  the  axial  component  of  vector  potential  for  the  azimuthal 
self-magnetic  field,  and  “a*  8s’  and  yg  are  constants. 

The  density  profile  of  the  solid  beam  described  by  Eq.  (6) 
is  given  by^ 


n°(r) 


0  <  r  <  R 


s 


(7) 


otherwise 


where  the  radius  of  the  solid  beam  is  determined  from  R  = 

s 

2  +  - 

v/(ui  -  w  ) (w  -  w  ),  the  effective  thermal  velocity  v„  is  defined 
s  s  s  s  s  » 


by  v 


2c2(ys  -  Ys)/Ys 


2  cs 


1  ± 


2  a)2 
ps  \ 

2  2  I 
Y  u 

s  cs 


1/2 


(8) 


are  the  laminar  rotational  frequencies  in  the  fast  (+)  and  slow  (-) 

2  -1/2 

rotational  equilibria,  e  *  sgn(e  ),  y  =  (1  -  3  )  ,  and  w  * 

s  s  s  s  cs 


eZ  B„/y  m  c  and  <L 
s  0  s  s  ps 


2  1/2 

(4irn  e  /y  m  )  are  the  cyclotron  and  plasma 

S  5  S  S 


frequencies,  respectively,  of  the  solid  beam  particles.  Obviously, 


the  equilibrium  solution  of  the  solid  beam  exist  only  for  the  rotational 
frequency  a)g  satisfying  wg  <  u>g  <  ug. 
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For  present  purposes,  we  also  specialize  to  the  sharp-boundary 
density  profile  of  the  hollow  electron  beam  given  by 


"h(r)  '  | 


n^  *  const.,  <  r  <  R2  , 


otherwise 


Since  the  mean  velocity  of  an  electron  fluid  element  in  the  hollow 
electron  beam  is  specified  by  ^(r)  ■  Vh0^r^8  +  ®hc.|z '  the  8tea^y~ 
state  equation  of  motion  in  the  azimuthal  direction  can  be 
expressed  as 


VVh9(r)/r  "  e[Er(r)  ‘  ShB8(r)  +  Vh8(r)B0/c]  *  (10) 

from  Eq.  (4).  In  Eq.  (10),  E^(r)  and  B^(r)  are  the  equilibrium  radial 

r  9 

electric  and  azimuthal  magnetic  fields,  and  e„  and  &  are  unit  vectors 

"  XQ  xz 

in  the  radial  and  azimuthal  directions.  Making  use  of  the  density 

profiles  in  Eqs.  (7)  and  (9),  it  can  be  found  that  the  field  combination 

E°(r)  -  S.  B^(r)  in  Eq.  (10)  is  expressed  as 
r  n  0 


Er(r)  -  6hB0(r)  “  *  17  U0(r)  "  8hAz(r)1 


2ire  n  (1  -  B  0.  )  — - 

ss  s  n  r 


27reAh 


~  R£/r)  , 


for  R^  <  r  <  R^.  Substituting  Eq.  (11)  into  Eq.  (10),  and  making  use 
of  Eqs.  (2)  and  (3),  we  obtain  the  rotational  frequency  w^(r)  of  the 


hollow  beam 


“h(r>  ■  “d(1  ■  7  *  y 


7  Z  Rc1 

h  (1  ”  ssBh)7  • 


where  the  diocotron  frequency  is  defined  by 


.  2  2 

u.  *  li  ,/  2y ,  j  , 
D  ph  h  ch 
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STABILITY  ANALYSIS  FOR  COUPLED  TRANSVERSE  PERTURBATION 

In  this  section,  we  make  use  of  the  linearized  Vlasov- fluid  and 
Maxwell  equations  to  obtain  the  dispersion  relation  for  coupled 
transverse  oscillation  of  an  intense  charge  beam  propagating  through  a 
relativistic  hollow  electron  beam.  We  adopt  a  normal-mode  approach 
in  which  all  perturbations  are  assumed  to  vary  with  time  and  space 
according  to 

6*(£.t)  ■  $(r)exp{i(Jte+kz-ti)t)}  , 

where  u  is  the  complex  eigenf requency  with  Imo>  >  0,  1  is  the  azimuthal 
harmonic  number,  and  k  is  the  axial  wavenumber.  The  present 
stability  analysis  is  carried  out  for  long  axial  wavelength 
perturbations  satisfying 

kV  «  a2  +  i)  ,  (w) 

c 

where  Rc  is  the  radius  of  the  conducting  wall.  Moreover,  the 
stability  properties  are  investigated  for  low-frequency  perturbation 
satisfying 

jwRc|2/c2  «  a2  +  1)  -  (15) 

Within  the  context  of  Eqs.  (14)  and  (15),  it  is  straightforward 

to  show  that  the  Maxwell  equations  for  perturbed  fields  can  be 
9 

approximated  by 

(?!? 't7-7)*(rt  •  <16) 

and 
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/  1  3  3  sr  \  ,  4ir  -T/  . 

(■?37r  3T“"2;A(r)  “  -r J(r)  * 


(17) 


where  $(r)  and  p(r)  are  the  perturbed  electric  potential  and  charge 
density,  respectively,  and  A(r)  and  J(r)  are  the  axial  components  of 
the  perturbed  vector  potential  and  current  density,  respectively. 
Components  of  the  perturbed  fields  can  be  expressed  in  terms  of 

A  A 

<j>(r)  and  A(r)  as 


E  (r)  -  -  (3/3rH(r)  ,  E.(r) 

r  o 


(18-/r)d> (r)  , 


B  (r)  «  -  (3/3r)A(r)  ,  B  (r)  -  (ii./r)A(r)  , 
“  r 


(18) 


E  (r)  -  -i[k$(r)  -  (io/c)A(r)  ]  . 
z 

The  perturbed  charge  ph(r)  and  current  Jh(r)  densities  contributed 
by  the  hollow  electron  beam  can  be  calculated  by  linearizing  the 
macroscopic  fluid  descriptions  in  Eqs.  (4)  and  (5).  Defining 
the  effective  perturbed  hollow  potential 

^h(r)  «  #(r)  -  BhA(r)  ,  (19) 

and  carrying  out  a  straightforward  algebraic  manipulation,  we  obtain 
the  hollow  beam  portion  of  the  perturbed  charge  density^ 


ph(r)  "  Jh(r)/ShC 


!>|ih(r) 

4irr 


(20) 


r  ,  _  .  ,  .  ,  3r  ph 

wch[ai  -  kfi^c  -  (r)  ] 


V 


2  2  0 

where  u  .  (r)  *  4ire  n,  (r)/y.m  is  the  plasma  frequency-squared  of  the 
pn  n  n 

2 

hollow  beam,  and  its  derivative  (3/3r)w  ,  (r)  can  be  expressed  as 

ph 
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fc“ph(r)  ‘  >t6(r  "  Rl)  '  5(r  -  R2)]  * 


from  Eq.  (9).  On  the  other  hand,  making  use  of  the  linearized  Vlasov 

equation,  the  perturbed  charge  density  pg(r)  contributed  by  the  9olid 
9 

beam  is  given  by 


ps<r>  -  Vr)/esc 


*4  (r) 

V  • 


where 


i/*g(r)  *  $(r)  -  SgA(r)  , 


is  the  effective  perturbed  solid  potential,  and 

l  ^TcfWr 

iv  \u  -  u  /  n*0 

s  k  s  s 

x  M  ~  *Ms  “  k6gc  /US  -  MSU 

in  -  k8  c  -  Hu  -  n(u*  -  u  )  \ui  -  m*/  j 


oj  -  k8  c  -  Zai 
s  s 


is  the  effective  susceptibility  of  the  solid  beam. 

Substituting  Eqs.  (20)  and  (21)  into  Eqs.  (16)  and  (17),  and 
making  use  of  the  definitions  of  the  effective  potentials 

*  a 

ifi  (r)  and  \ p  (r)  give  the  coupled  eigenvalue  equations 
n  s 


\  r  3r 


3r  r  3r 


')*b(t) 


i?>h(r) 


H(r)  [6 (r  -  R^  -  6(r  -  R^)] 


i<Pa(r) 


O.  -  SgSh)S6(r  -  Rg)  , 


(r^fT-MV^ 


(r) 


S6(r  -  Rg) 
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**h(r) 


(1  -  Ss6h)H(r)[6(r  -  R^  -  6  (r  -  R£)]  , 


(25) 


where  the  effective  susceptibility  of  the  hollow  beam  is  defined  by 


H(r)  -  “ph^ch^  “  ^hc  ~  iwh(r)  1  •  (26) 

It  is  evident  that  the  right-hand  sides  of  Eqs.  (24)  and  (25) 
are  equal  to  zero  except  at  the  surface  of  the  beam  boundaries 
(r  ■  Rg,  r  »  R^,  and  r  =*  R2) .  Therefore,  the  eigenvalue  equations 
(24)  and  (25)  can  be  reduced  to  r-1(3/3r) (r3^/3r)  -  (Z2/r2)i}ij  »  0, 
except  at  r  ■  Rg,  r  *  R^,  and  r  *  R2*  The  eigenfunction  satisfying 

Eqs.  (24)  and  (25)  are  expressed  as 

J. 


-l 


Vr>  *  + 

Djr  +  E  r  , 

Fj  (r*  -  R**r~*), 


0  <  r  <  Rg  , 
Rg  <  r  <  Rx  , 

Ri  <  r  <  R2  ’ 

<  r  <  R„ 


(27) 


R2  '  *■  '  *'c 

where  the  subscript  j  »  h  and  s  denote  the  hollow  and  solid  potentials, 

and  the  constants  Aj  -  Fj  can  be  related  by  integrating  Eq.  (27) 

across  the  discontinuities  at  r  *  R  ,  R  ,  and  R_.  After  carrying  out  a 

si 

straightforward  algebraic  manipulation,  we  obtain  the  dispersion 
relation  for  coupled  transverse  oscillation 


rh(«,k)  •  r  (w,k) 


Vsa  -  M 


H(R, 


Vv  V 

:(-S)(-f)l 


s°h' 

21, 


H(R2)  Rx 


21 


.21 


M 


1  - 


.21 


.21 


! 1  -  “77  )  h  (28) 


where  the  dielectric  functions  of  the  hollow  and  solid  beams,  T^(w,k)  and 

T  (w,k),  are  defined  by 
s 
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rh(u,k) 


(29) 


and 


r  (<u,k)  -  —^2  (i  - 

2y!  \  R‘ 


2*> 

s 

2  Jt 
c  ' 


1  . 


(30) 


Equation  (28)  is  the  dispersion  relation  used  in  the  remainder  of 
this  paper,  and  can  be  used  to  investigate  stability  properties  for  a 
broad  range  of  system  parameters  of  experimental  interest. 

In  the  limit  of  S  ■*  0,  we  recover  the  dispersion  relation  of  the 
diocotron  instability^- 


rh(u,k)  -  0  ,  (31) 

for  a  relativistic  hollow  beam  from  Eq.  (28).  On  the  other  hand, 
in  the  absence  of  hollow  beam  (i.e.,  H  ■+■  0)  ,  Eq.  (28)  reduces  to 

r8(«,k)  -  0  ,  (32) 

which  is  identical  to  the  result  obtained  by  Uhm  and  Davidson^  for 

a  relativistic  nonneutral  electron  beam.  A  careful  examination 

of  the  right-hand  side  of  Eq.  (28)  shows  that  the  coupling  between 

the  hollow  and  solid  beam  modes  is  minimum  whenever  3g  *  8^*  Moreover, 

for  a  solid  beam  with  its  radius  Rg  satisfying  R  <<  R]_,  the  fundamental 
•  s 

mode  (£,  »  1)  perturbation  is  the  most  unstable  mode  in  the  coupled 
transverse  oscillations. 

An  important  parameter  that  characterizes  the  importance  of  kinetic 
effects  for  solid  beam  is  the  ratio  of  the  characteristic  thermal 

Larmor  radius  (rT  *  v  /u  )  to  the  beam  radius  R  .  Making  use  of  the 

LS  S  CS  S 
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beam  radius  definition  in  Eq.  (7),  it  is  straightforward  to  show  that 


rf  /R^  *=  (w  u+)  (id  -  u>  )/w^ 
Ls  s  s  s  s  s  cs 


where  u'  is  defined  in  Eq.  (8),  Evidently,  the  cold-fluid  limit, 

2  2 

where  the  beam  motion  is  approximately  laminar  with  r  <<  R  , 

LS  S 

+ 

corresponds  to  rotation  velocities  that  satisfy  ws  3  A  careful  exam¬ 


ination  of  the  expression  for  S  shows  that 
lim  +  S(u,k) 


18 


~  2 
u 

_£S_ 


s  s 


(w  -  kg  c  -  -  k6  c  -  iu-)  +  e  tu  +  2u>  ] 

s  s  s  s  s  cs  s 


(33) 


Although  the  finite  Larmor  radius  effects  of  the  solid  beam  particles 
18 

are  important  on  the  stability  properties,  the  numerical  analysis  of 

the  dispersion  relation  in  Eq.  (28)  is  restricted  to  the  cold  particles 

characterized  by  £  /R  -*•  0.  Moreover,  in  the  remainder  of  this  paper, 
l*s  s 

we  assume  that  the  solid  beam  particles  are  in  a  slow  rotational 
equilibrium,  i.e.,  w  =  a>  . 
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STABILITY  PROPERTIES  OF  A  SOLID  ELECTRON  BEAM 


In  this  section,  we  investigate  the  stability  properties  for 
coupled  transverse  oscillations  between  the  solid  electron  beam  and 
the  hollow  electron  beam.  Since  the  cyclotron  frequency  of  the  solid 
beam  electrons  is  much  higher  than  its  rotation  frequency 


ID 

S 


0) 

s 


.2  n  2 
ui  /  iy  u 
ps  s  cs 


(34) 


or  the  eigenfrequency  (i.e.,  w  >>  u  or  u>  >>  ai) ,  the  effective 

cs  s  cs 

susceptibility  S(ui,k)  of  the  solid  beam  in  Eq.  (33)  can  be  approximated  by 

.2 

ai 

S(ui,k)  - - £5 - z -  .  (35) 

(u)  -  kS  c  -  4uj  )u> 

s  s  cs 

Defining  the  normalized  Doppler-shifted  eigenfrequency  0  by 


fl  ■  (u  -  k8sc)/uiD  ,  (36) 

and  substituting  Eq.  (35)  into  Eq.  (28),  we  calculate  the  normalized 

growth  rate  0^  *  Imft  and  Doppler-shifted  real  oscillation  frequency 

»  ReD  numerically  from  Eq.  (28)  for  a  broad  range  of  system 

parameters,  8.,  8,  ,  B  ,  n  / n,  ,  R  /R  ,  R./R  ,  and  Ro/R  .  However,  we 
nssnscl.£  *  c 

emphasize  that  the  normalized  growth  rate  JL  calculated  in  this  section 

2  2 

is  independent  of  the  parameter  “ph/(1)ch- 

Shown  in  Fig.  1  is  plots  of  the  normalized  growth  rate  »  lmfl 

versus  kc/w^  for  R^/Rc  =  0.857,  “  0.939,  Rg/Rc  *  0.2,  8^  *  0.968  (y^“^)  * 

4  *  1,  n  /n,  *  0.1,  and  several  values  of  S=-  For  a  small  radius  solid 
s  h  s 

19 

beam  (R  /R  ■  0.2  <<  1)  consistent  with  the  present  experiment, 

&  C 

the  4*1  perturbation  is  the  most  unstable  mode.  In  these  physical 
parameters,  the  growth  rate  of  instability  for  a  coupled  transverse 
oscillation  is  a  substantial  fraction  of  the  diocotron  frequency 
thereby  severely  limiting  the  solid  beam  propagation.  The  axial  wavenumber 


21 


NSWC  TR  81-391 


and  the  real  oscillation  frequency  corresponding  to  the  maximum  growth 
rate  of  instability  satisfy  (kc/w^l  <  1  and  |n  |  *  0.1,  respectively. 
However,  we  also  note  from  Fig.  1  that  small  growth  rate  in  range  0.5  < 
kc/a)jj  <  2  is  a  residual  influence  of  the  familiar  diocotron  instability. 
Therefore,  we  conclude  that  the  growth  rate  of  the  diocotron 
instability  is  much  less  than  that  of  the  coupled  transverse  oscillation 
for  counter-streaming  beams. 

The  dependence  of  stability  properties  on  axial  velocity  of  the 

solid  beam  is  illustrated  in  Fig.  2,  where  the  normalized  maximum 

growth  rate  0®  is  plotted  versus  Bg  for  the  parameters  identical  to 

Fig.  1.  The  maximum  growth  rate  of  instability  almost  linearly  reduces 

to  zero  as  the  value  8  Increases  from  B  »  -1  to  B  *  B,  *  0.968. 

s  s  s  h 

In  this  regard,  it  is  evident  that  the  1*1  mode  perturbation  is 

stable  when  the  solid  and  hollow  beams  have  a  same  axial  velocity. 

2 

This  is  consistent  with  the  previous  study.  Of  considerable  interest 
for  experimental  application  is  the  stability  behavior  for  specified 
values  of  8^  and  8g.  Typical  result  is  shown  in  Fig.  3  where  the  maximum 
growth  rate  is  plotted  versus  ^g/n^  for  0g  *  -0.943  (yg  *  3)  and 
parameters  otherwise  identical  to  Fig.  1.-.  Obviously  from  Fig.  3,  the 
maximum  growth  rate  of  Instability  reduces  drastically  with  decreasing 
value  of  the  density  ratio  hg/n^.  We  therefore  conclude  that  a  . 
solid  beam  with  a  substantial  density  reduction  can  propagate  through  a 
counter-streaming  hollow  beam,  without  a  severe  growth  of  unstable 


transverse  oscillations. 
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ION  RESONANCE  INSTABILITY  IN  A  HOLLOW  ELECTRON  BEAM 

18 

The  Ion  resonance  instability  is  one  of  the  fundamental 
instabilities  that  characterize  a  relativistic  nonneutral  plasma 
system  with  both  ion  and  electron  components.  In  this  section,  we 
investigate  ion  resonance  stability  properties  of  a  solid  ion  beam 
propagating  through  a  hollow  relativistic  electron  beam.  The  stability 
analysis  in  this  section  is  restricted  to  Rg  <  which  requires  a  low 
density  ion  beam  satisfying 


2  2 
y  u> 
s  cs 


(37) 


For  the  ion  resonance  instability  in  a  high  density  ion  beam  with 

2  2  2  9 

2w  >  y  lit  ,  we  urge  the  reader  to  review  the  previous  study, 
ps  s  cs 

Defining  the  normalized  Doppler-shifted  eigenf requency  !1  by  fl  »  (w  -  kS  c)/<d 

5  U 

in  Eq.  (36),  we  have  obtained  the  growth  rate  and  real  oscillation 


frequency  numerically  from  Eq.  (28)  for  a  broad  range  of  system 
2  2 

parameters  <L  ,  /w  .  ,  R  /R  ,  8  ,  and  n  /n,  .  In  this  section, 
ph  ch  s  c  s  s  h 

we  summarize  the  essential  features  of  the  stability  studies.  The 


numerical  analysis  is  restricted  to  the  choice  of  parameters 

B.  -  0.995  (y.  -  10),  £2./<A  -  0.1,  R./R  -  0.857,  R,/R  -  0.939  and 
h  1  h  ph  ch  i  c  z  c 

mg/m  *  1836  (proton  beam). 

The  dependence  of  stability  properties  on  axial  wavenumber  k  is 
illustrated  *in  Fig.  4,  where  (a)  the  normalized  growth  rate  and  (b) 
Doppler-shifted  real  oscillation  frequency  are  plotted  versus  kc/w^ 
for  Rs/Rc  *  0.8,  Bg  *  0.1,  ng/hh  •  0.01,  and  the  azimuthal  harmonic 
number  4  ■  1  and  2.  In  Fig.  4(b),  the  real  oscillation  freauency  is 
plotted  only  for  the  ranges  of  kc/u^  corresponding  to  instability  (lm  H  >0). 
Several  features  are  noteworthy  in  Fig.  4.  First,  the  maximum  growth 
rate  of  instability  is  many  times  of  the  diocotron  frequency,  thereby 
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severely  limiting  the  ion  beam  propagation  through  a  relativistic  hollow 
electron  beam.  Second,  as  shown  in  Fig.  4(a),  where  R  /R  2  1, 

S  X 

perturbations  with  azimuthal  harmonic  number  4  >  2  are  the  most  unstable 

mode.  Finally,  we  note  from  Fig.  4(a)  that  the  growth  rate  of  the 

dlocotron  instability  in  the  range  kc/w^  >  10  is  much  less  than  that 

of  the  ion  resonance  instability.  Shown  in  Fig.  5  are  plots  of  the 

normalized  growth  rate  versus  kc/w^  for  4-1,  several  values  of 

Sg  and  parameters  otherwise  identical  to  Fig.  4.  As  expected,  the 

growth  rate  of  instability  is  a  decreasing  function  of  the  parameter  B  . 

s 

In  order  to  illustrate  the  dependence  of  stability  properties  on 
the  solid  beam  radius,  we  present  the  normalized  growth  rate  versus 
n  /fl,  in  Fig.  6  for  several  values  of  4,  and  (a)  R  /R  -  0.8 

S  il  sc 

and  (b)  R  /R  *  0.4.  It  is  evident  from  Fig.  6  that  the  instability 
growth  rate  increases  significantly  as  the  solid  beam  radius  Rg  is 
increased.  Moreover,  it  is  also  noted  from  Fig.  6(b)  that  for  small 
beam  radius,  the  4=1  mode  is  the  most  unstable  mode.  For  example, 
the  growth  rate  ^of  the  4  =  4  perturbation  in  Fig.  6(b)  is  that  of  the 
typical  diocotron  instability.  After  a  careful  examination  of  Figs. 

6(a)  and  (b),  we  finally  note  that  for  a  specified  value  of  (n  /h.  )(R  /R  ) 
corresponding  to  the  total  ion  current,  the  growth  rate  of  the  coupled 
transverse  oscillation  reduces  considerably  as  the  ion  beam  radius  is 


decreased. 
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CONCLUSIONS 


In  this  paper,  we  have  investigated  the  stability  properties  of 
coupled  transverse  oscillations  of  a  solid  particle  beam  propagating 
through  a  relativistic  hollow  electron  beam.  The  analysis  was  carried 
out  within  the  framework  of  a  hybrid  (Vlasov-fluid)  model  in  which  the 
hollow  beam  electrons  are  described  as  a  macroscopic  cold  fluid  and  the 
solid  beam  particles  are  described  by  Vlasov  equation.  Equilibrium 
properties  and  the  basic  assumptions  were  presented  in  Sec.  II. 

Stability  analysis  of  coupled  transverse  oscillations  was 

carried  out  in  Sec.  Ill  for  long  axial  wavelength  and  low-frequency 

perturbations.  A  closed  algebraic  dispersion  relation  in  Eq.  (28) 

was  obtained  from  the  solid  and  hollow  beams  with  sharp-boundary 

density  profiles  in  Eqs.  (7)  and  (9).  One  of  the  most  important 

features  in  the  analysis  is  that  the  coupling  between  the  hollow  and 

solid  beam  modes  is  minimum  whenever  both  beams  have  a  same  axial  velocity. 

Stability  properties  for  coupled  transverse  oscillations  between 
the  solid  electron  beam  and  the  hollow  electron  beam  was  investigated 
in  Sec.  IV.  We  found  that  the  typical  growth  rate  of  the  transverse 
oscillation  is  a  substantial  fraction  of  the  dlocotron  frequency 
thereby  severely  limiting  the  solid  beam  propagation  through  the  hollow 
electron  beam.  The  growth  rate  of  the  conventional  diocotron  instability 
is  much  less  than  that  of  the  coupled  transverse  oscillations  for  counter¬ 
streaming  beams.  Finally,  the  ion  resonance  stability  properties  of  a 
solid  ion  beam  was  investigated  in  Sec.  V  for  a  low  density  case.  It 
was  found  that  the  maximum  growth  rate  of  the  ion  resonance 
instability  is  many  times  of  the  diocotron  frequency.  Moreover,  the 
instability  growth  rate  increases  significantly  as  the  solid  beam  radius 
Rg  is  increased. 
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£=\.  £h  =  0.968,  ns/nh  =  0.1,  Rs/Rc  =  0.2 
R,/Rc  =  0.857,  R2/Rc  =  0.939 


/3S=0 


@s  =  0-5 


A=-0.5 


kc/oj, 


Fig.  1  Plot  of  normalized  growth  rate  £1^  versus  kc/w^  for  electron- 

electron  interaction,  R^/Rc  *  0.857,  *  0-939,  Rg/Rc  ■  0.2, 

8.  ■  0.968,  4  ■  1,  n  /n.  ■  0.1,  and  several  values  of  6 
n  s  n  s 
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Fig.  3  Plot  of  normalized  maximum  growth  rate  ft™  versus  ng/nh  for  3s  --0.943 
<YS  *  3)  and  parameters  otherwise  identical  to  Fig  1 
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l-l  and  2 
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kc/ uj  q 


Fig.  5  Plot  of  normalized  growth  rate  versus  kc/ui^  for  l  •  1, 
several  values  of  Bs  and  parameters  otherwise  identical  to 
Fig.  4 
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